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We discuss the general framework for a perfect continuum medium in cosmology and show that 
an interesting generalization of the fluids normally used is for the medium to have rigidity and, 
hence, be analogous to an elastic solid. Such models can provide perfect, adiabatic fluids which are 
stable even when the pressure is negative, if the rigidity is sufficiently large, making them natural 
candidates to describe the dark energy. In fact, if the medium is adiabatic and isotropic, they 
provide the most general description of linearized perturbations. We derive the equations of motion 
and wave propagation speeds in the isotropic case. We point out that anisotropic models can also 
be incorporated within the same formalism and that they are classified by the standard Bravais 
Lattices. We identify the adiabatic and isocurvature modes allowed in both the scalar and vector 
sectors and discuss the predictions they make for CMB and matter power spectra. We comment 
on the relationship between these models and other fluid-based approaches to dark energy, and 
discuss a possible microphysical manifestation of this class of models as a continuum description of 
defect-dominated scenarios. 
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I. INTRODUCTION 



> 

£-■*. , There are four stress-energy components in the presently preferred standard cosmological model [3, [E El : the 
baryons which constitute the visible universe, radiation (photons and neutrinos), cold dark matter (CDM) and some 
unknown component, often called dark energy, which gives rise to the cosmic acceleration [H, H, IE 0, H|. Within the 
codes [E [T3| used to make predictions for specific models, the CDM is modelled as a pressureless, perfect fluid and 
the radiation components are modelled as black-body radiation gases. However, many of the qualitative features of 
• the observed power spectra, for example, the acoustic peaks in the angular power spectrum of the cosmic microwave 

CLf background (CMB), can be derived by ignoring the higher order moments of the photon and neutrino distributions 
and treating them as a perfect fluid with P = p/3, where P is the pressure and p is the density. Hence, most of the 

$H ' important components of the universe are well approximated by perfect fluids. 

The microphysical origin of the cosmic acceleration is still a mystery. The simplest, and probably most popular, 
explanation is a cosmological constant, although this has some well documented fine-tuning problems |Tll.[l2l|. Various 
ideas exist, collectively known as dark energy (for a recent review of dark energy models, see ref. |l3|), whereby the 
additional stress-energy component is provided by a slowly rolling scalar field 0, EE EH , known as as Quintessence, 
or a lattice of topological defects [IE EE EE H3|- Alternative explanations require the modification of gravity by 
the inclusion of non-minimal coupling between matter and gravity [2lj . extra-dimensions (for example, ref. [22| ) and 
modifications to the Einstein-Hilbert action [IE [24], [H| . 

In this paper we shall consider a generalization of the perfect continuum fluid approach to describe dark energy. We 
have already pointed out that perfect fluids can be used to describe the essential properties of the radiation and CDM 
in the universe, and it seems reasonable to consider the possibility that the dark energy can also be understood in 
the same way. This will require us to go back to the fundamentals of how formulate a generalized continuum medium 
in General Relativity and derive equations for the perturbations in the medium. These are essential in computing 
accurate predictions for observed power spectra [IE HE HE H3 • We will show that the most obvious generalization of 
a perfect fluid is to allow for rigidity of the medium in a way analogous to a continuum elastic solid [19( , and that this 
can be stable if the rigidity is sufficiently large even if the pressure is negative. For most cosmological observations 
(for example, the CMB or large-scale structure measurements) we only need to consider the linearized regime. We 
will argue that the generalization of fluids to include rigidity, is the most general possibility for an adiabatic medium 
at linearized order. 

In many ways the concept of dark energy is similar to the idea of the Aether postulated in the late 19th century: 
there is something about the laws of physics which appears awry and we postulate a medium to try and solve the 
problem (65|. The approach we propose is to formulate generalized properties of such a medium, define ways of 
computing power spectra and then compare them with observations (see ref. [28j for the status of this kind of model 
after the first year WMAP data). 
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The idea of generalized dark matter /energy has been considered previously by a number of authors [lj| [2(| [26|, 
I271 . [29l [30l . [3lj . We will attempt to discuss how the approaches suggested by others relate to those presented here. 
The broad difference between our approach and those taken previously is that we will derive the equations describing 
perturbations in the medium from a set of well defined physical assumptions. These equations are closed and there is 
no freedom, except the strength of the rigidity, to play with. 

The original derivation of the equations presented here , was motivated by the idea that the dark energy could 
be a lattice of topological defects (cosmic strings or domain walls) formed at a low energy phase transition. We will 
explain how such a lattice can provide a possible microphysical model for an elastic dark energy model. However, we 
believe that the formulation of the elastic dark energy models is more general and should not be thought of as being 
necessarily linked to these specific type of models which predict very specific values for w — P/ p. In particular, the 
elastic dark energy models include CDM and a cosmological constant as limiting cases allowing them to provide an 
interesting phenomenology. 

The paper is organized as follows. In section |TT] we discuss the formulation of a generalized medium and derive 
the equations of motion. We then go on to identify the various adiabatic and isocurvature modes in section [TTT] and 
present the cosmological signatures expected in the CMB and matter power spectra in section IIVI In section [V] we 
explain the basic ideas of defect-dominated scenarios for dark energy and explain how they relate to the elastic dark 
energy models. 



II. COSMOLOGICAL DYNAMICS OF A GENERALIZED MEDIUM 



In Section UTAl we review the medium representation concept and the scheme for specifying perturbations in a 
general medium. In III Bl we apply this formalism to compute the cosmological equations of motion for an isotropic 
elastic medium and in III d evaluate the propagation speeds of perturbations. We then decompose the perturbations in 
terms of harmonic basis functions in III Dl and III El to give the full set of Einstein and energy- momentum conservation 
equations. In Section Hi Fl wc compare our results to other treatments of generalized fluids, and in III Gl discuss briefly 
how the formalism can be applied to anisotropic perturbations. 



A. Medium Representation Concept 

This section provides an overview of work on relativistic elastic media in the context of neutrons stars by Carter 
and others [32L I33L l34l . [35| , and details how this elegant formalism to describe the mechanics of a generalized medium 
can be applied to cosmology. In Newtonian theory it is natural to define the properties of a medium, such as the 
density and pressure, in terms of coordinates of a three-dimensional Euclidean space at some given instant of time. In 
General Relativity, however, there is no general time-slicing of the four-dimensional space-time manifold M. which can 
be used to specify the material state. Therefore, it is necessary to consider the projection V: A4 — > Ti of A4 onto a 
three-dimensional manifold Ti whose elements represent particles of the medium [32]]. The inverse image V~ 1 {X) C M 
of a point X G Ti can then interpreted as the worldline of the particle represented by X. The inverse image determines 
a one-to-one mapping between medium tensors defined on Ti and the set of tensors on M. which are orthogonal to the 
congruence of wordlines. This mapping is important as it allows the intrinsic properties of the material medium to 
be defined on Ti, whilst the spacetime evolution is described by tensor fields on M. Moreover, it allows the definition 
of spacetime tensors in terms the tensors defined on Ti. We will use Greek indices p, v... to label tensors on M and 
Roman indices A. B... to label tensors on Ti. 

We define jab to be the metric on Ti which quantifies the strain of the medium. At each point of Ti, we will assume 
that there are functions determining the density p and pressure P AB which can be expressed in terms of jab (one 
of the properties of a perfect medium) and, most probably, related by an equation of state P AB (p). Under these 
assumptions one can show that 

where (7) is the determinant of 7ab, and by taking a second derivative of p with respect to 7ab, one can define 

E ABCD = _ 2| pl/2 3_ (| | l /2pCD) (2) 
0(1AB) 

the classical elasticity tensor which satisfies E ABCD = E( AB K CD ) = E CDAB . We will see that defining these two 
tensors is sufficient to understand linearized perturbations of the medium. 



3 




FIG. 1: Evolution of a material along the flow lines u M . Physical quantities are defined on the material manifold TL in terms 
of internal co-ordinates A,B..., and there is a direct mapping of spacetime tensors orthogonal to the flow and tensors on TL. 
Parallel components of the spacetime tensors do not contribute to the material projection. 



In order to make correspondence with the spacetime A4, we define to be a bi-tensorial projection operator 
which projects a tensor in TL into A4, or vice- versa, then one can define the projected metric tensor 7 M „ by 

7m" = VffilAB ■ (3) 

The field of flow vectors u M tangent to wordlines in M. are normalized by the condition 

«\ = -l, (4) 

and the projection is orthogonal to them, that is, u^V^ — 0. This allows us to construct 7 M „ as an orthogonal 
projection tensor 

7/U> = + U P- U v ' ( 5 ) 

providing an alternative, but equivalent, definition for j^. It acts as the tensor which projects other tensors onto the 
tangent subspace orthogonal to the wordlines. For example, any vector v p can then be decomposed into components 
orthogonal and parallel to the flow by =±v p ' + v^u^, where ju M = J^v" and = —u^v^. This allows physically 
relevant spacetime tensors to be decomposed into the parallel component, which does not contribute to the material 
projection, and the orthogonal part, which does. This framework is summarized schematically in Fig.[TJ 
One can use the projection tensor to construct the spacetime pressure tensor, P^ v , from that on TL 

P»u = VffiPAB , (6) 
which satisfies u^P^ — 0. The same can be done for the elasticity tensor E^ vp<T which satisfies 

JTj^vpa _ j^{pv){pa) _ j^pcrp,v E^ 9 " U a = (7) 

One must also specify how perturbations in TL are related to those in M. in order to deal with dynamics. The 
convected differential d[ ] is an important quantity, as it is the spacetime tensor corresponding to the Lagrangian 
(wordline preserving) material variation on TL. In the case of a vector, there is the bijection d[t>' 1 ] < — > {<5l aP A , <5lv"}, 
where <5l is the Lagrangian variation. The convected differential can be evaluated in terms of Lagrangian differentials 
on M. by projecting the orthogonal part of the perturbation. This gives the explicit relationship [34[ 

d[T£;;] = * L T£" + T^;;u"^u v ~ T£;;;u»S L u p , (8) 

where T£"' is a general mixed tensor. 

The next stage is to relate the material variations in terms of spacetime tensors. The most direct way to do this 
is to use the convected derivative, which relates Lagrangian and material variations. Applying ([5]) to (P) and ((5J) one 
obtains [H| 

hp = ~ (P^ + fr/n^d^, (9) 
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J*" = -^(E pl/pa + P pv Y° - ±P p{p u v) u a )5^g p(7 , (10) 



where extra terms are induced in (|T0|) due to ([8]). 

We can now begin to discuss the dynamics of an elastic medium. A perfect elastic medium is defined by the condition 
that the energy-momentum tensor T pv is a material function of the metric tensor with respect to the flow field [35|, 
meaning that it takes the form 

T pv = pu^u v + P pv . (11) 

The variation in the energy-momentum tensor can be obtained by using (jlip in conjunction with ([9]) and (|10|) to 

give 

5 L T^ = -l(W^<>° +T^g?°)5 h g pa , (12) 
where the non-orthogonal elasticity tensor W pvpa can be decomposed as [36j j 

w ^pa = E ^pa + p^ u p u a + ppttftf _ p^ u ^ _ P^ U V U P - P Vp U P U° - P V ° U P U p - pU P U U U P U a , (13) 

and has the same symmetry properties as the ordinary elasticity tensor. This allows the variations to be conveniently 
written in terms of functional derivatives of the Lagrangian with respect to the metric via 

T^ = -2\g\-^^—{\g\^C), (14) 
W^ p ° = 4| 3 r 1 /2-^-^(| 5 |i/2 £) = -2\g\-^-^{\g\^Tn , (15) 

Oh9pa OLg^u Ol^gpa 

where \g\ is the determinant of the metric. 

It is often more convenient to describe perturbations being fixed with respect to some background space. If the 
vector field ^ p is the infinitesimal displacement of the wordlines with respect to the fixed background space, then the 
difference between Lagrangian <5l and fixed (Eulerian) variations <5e is given by 

5 L = S E + £ € , (16) 

where is the Lie derivative. Such a displacement could be removed by using a mapping of the perturbed space 
onto the background space, but it will be more convenient to set the gauge by some other means. In the following 
section, for example, we use the synchronous gauge to derive the cosmological equations of motion. For the particular 
case of the metric tensor then the relationship (Til)]) gives 

<Sl 9pu = Se 9pu + 2V( Jll &,) , (17) 

where <5e 9pv is the Eulerian variation of the metric tensor. 

The equation of motion for the vector field £ M gives the complete system of equations for the perturbations. 
Evaluation of the Lagrangian variation <5l(7 a1 1/ V /j T' a "') gives 

{A* v f) - ( P1 p P + P%)u v u°)5^ va + lr p l\l P ASL g a p)V T E p ™° = (18) 

(P^if _ Ip-^M _ 2A«% T )vP T u a + {pi p p + P» p )u p u v u°)8^g va , 

where dots denote covariant differentiation with respect to the flow (that is, u^V p ) and 

A pv * _ Eii v„ _ y^pw 5 ( 19 ) 

is the relativistic Hadamard elasticity tensor which forms the characteristic equation needed to evaluate the sound 
speeds in the medium (3?J ■ This tensor obeys the symmetry and orthogonality conditions 

A pv P a = A p* P ^ A pvpa u a = . (20) 

Eqn. (fl8|) takes the form of a wave equation for the displacement vector £ M . Since all components are orthogonal the the 
flow then the additional degree of freedom in £ M can be removed by imposing the orthogonality requirement = 0. 
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The remaining quantities required to evaluate (|18[) are the Lagrangian variation of the connection coefficients, given 
by 

Sir» va = V (l A g» a) - l - g va , (21) 

and the flow gradient tensor given by 

Vfu, = V, y u M + u^Uv . (22) 

B. Cosmological Equations of Motion for an Isotropic Medium 

The perturbed space-time metric takes the form 

9llv = a 2 {r) [r,^ + h^] , h» v = rf»jf v h afs , (23) 

and so the Eulerian component of the metric perturbation is given by Se g^v — or h^,. We make use of the synchronous 
gauge conditions [hoo — hoi — 0) to remove the remaining degree of gauge freedom in the Einstein field equations. At 
zeroth order the flow vector = a _1 (l,0, 0,0), and so the non-zero components of the displacement vector £ M are 
confined to the spatial part by the orthogonality requirement = 0. 

In the synchronous gauge the non-zero components of the Lagrangian variation of the connection are given by 

fc,r*oo = £ + ft£ , S L T°oi = nii, ,v , -lHi) X, I Uh, s ■ h,,. (24) 

SiFoi = d 3 C + \h) , fer^fc = djdkC - *,,uc + d u h\ } - ~d% k , 

where dots are now understood to denote derivatives with respect to the conformal time and H is the conformal time 
Hubble parameter. We can now obtain the equations of motion and perturbed energy-momentum sources for the 
cosmological fluids by inserting the appropriate pressure and elasticity tensor expressions for each component into 
(H2D and (HHD- 

1. Isotropic Perfect Fluid 

In an isotropic perfect fluid the pressure tensor is isotropic and is given in terms of the pressure scalar P by 

P" v = P^ v , (25) 
while the elasticity tensor is given in terms of the bulk modulus (3 by [12] 

E ^ P a = P ) 7 ^ 7 P<x + 2 p 7 MpyT)^ ( ( 26 ) 

and the bulk modulus is defined by 

f3=( P + P)^f. (27) 
dp 

Substituting these expressions into (|12[) . we obtain the components of the perturbed energy- momentum tensor 

SeT° q = (p + P) Ue + ^hj , (28a) 
SeT* = -(p + P)C, (28b) 
Sb'Tj = -0j U^ + ^h) , (28c) 

where h is the trace of the metric perturbation. From now onward we do not make the explicit distinction between 
Eulerian (fixed) and Lagrangian perturbations, as all relevant quantities are evaluated with respect to the fixed 
coordinates. Substitution of ([25]) and (|26|) into (fT8|) gives the equation of motion for the displacement vector as 

(P + P){€ + W) - 3pH? - f3{d l d ] e + d l h/2) = . (29) 
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2. Isotropic Perfect Elastic Medium 

We have pointed out that perturbations in an elastic medium can be specified by the pressure tensor P^ v and 
the elasticity tensor E pvpa ' . In the case of isotropy the pressure tensor is also given by ([251 since 7^ is the only 
isotropic tensor of rank two, upto a scaling. The case of rank four, required to describe the elasticity tensor, is more 
complicated; it is given by 

W vpa = A'f'f + 5 7 ^7°> , (30) 

where A and B are arbitrary parameters. In order to fit in with the definitions used in the previous section, we will 
define an additional shear contribution by 

E »upa = ^ypa + _ p^jP" + 2P 1 p(P l' 7> , (31) 

where the shear tensor obeys the symmetry and orthogonality conditions 

Y^p,vpa _ ^(nv)(pa) _ Yfopv YP ,Vp<T U a = (32) 

In an isotropic elastic fluid the shear tensor in terms of a single shear moduli p by 

YP>p° _ 2 p M" 7 ')" - ^7^7^ • (33) 

The contribution of the shear tensor to the elasticity tensor is zero in the perfect fluid case (p, — 0). We again 
substitute these expressions into (|12|) to obtain the components of the perturbed energy-momentum tensor 

ST° Q = (p + P) (diC + h) , (34a) 
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5T\ = -(p + P)C, (34b) 
ST'j = -tfjtf - (d k e + \h)j - /i(2%f ) + h)) . (34c) 

The equation of motion for the displacement vector £* is then given by 

[p + P)(e + HC) - 3/Wf - ft&dj? + d l h/2) - ni&djC + 3 + d J h) - d l h/3) = . (35) 

where (|34al I34b[ I34cp and ([33)) agree with the equations for an elastic medium defined in ref . [lj| . If the medium 
forms at some finite time then boundary conditions imply that h{j — > — hfj , where / refers to quantities defined at 
the time of formation of the medium. Since only the second-order variation of the Lagrangian is relevant for linearized 
perturbations we have shown that the pressure tensor, P pv , and elasticity tensor, E p,vpa ', specify the most general 
parameterization of perturbations in T pv under the assumptions discussed earlier. 



C. Evaluation of Sound Speeds 



This section provides a brief summary of how to compute the propagation speed and polarization directions of 
sound waves (small perturbations) in a perfectly elastic medium, which was originally developed in ref. [37J ■ A sound 
wave front is defined as the hypersurface across which the acceleration vector u^V^ii^ has a jump discontinuity (that 
is, at some coordinate point Xq the acceleration is a well defined function in the limits xt and x+ , but is not defined 
at Xq). The flow vector u p , the functions of state p, P Miy , E< J,vpa , the space-time tensor g^ v and the projection tensor 
7^1/ are continuous across this hypersurface, but the acceleration induces discontinuities in the first derivatives of p, 
P pu , E pvpa and 7 M!/ . 

At the wave front hypersurface u u V u u p — al^, where a is the amplitude of the sound wave and is the polarization 
vector of the wave front, which satisfies 

1% = 1, 1^=0. (36) 

If one introduces a propagation direction vector v p satisfying the same orthonormality conditions as the polarization 
vector, then a characteristic equation 



[v*0rT + P IM ')-QHl» = 0, 



(37) 
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can be constructed whose eigenvalues give the squared propagation speeds v 2 in the direction specified by and the 
eigenvectors are the polarization direction(s). The relativistic Fresnel tensor is defined by 

Q »u = A w VpVtr , ( 38 ) 

where the Hadamard tensor A^^ 17 is defined by (fl"9]) . The Fresnel tensor satisfies the symmetry and orthogonality 
conditions 

Q» V = Q^ V \ Q^u u = Q. (39) 



1. Isotropic Perfect Fluid 



The Fresnel tensor can be obtained by substituting the expression for the pressure tensor ([25]) and elasticity tensor 
(US]) into (HU) and then using J3BJ| to give 

= 0v"v u , (40) 

so the eigenvalue equation becomes 

[v 2 {p + PW J -^v v ]l v = 0. (41) 

This equation has a single solution in which the propagation direction is parallel to the polarization vector (l^ = v^) 
and is given by 

v 2 = c 2 = -J*— = ^ , (42) 
s p + P dp' v ' 

which corresponds to the longitudinal (scalar) sound speed. 



2. Isotropic Perfect Elastic Medium 



In the case of an isotropic elastic medium the Fresnel tensor is given by 



P + ''''<"'- //:■'"' 



so the eigenvalue equation becomes 



v 2 (p + p)Y u - mt^ - [P + iH v ^ vV 



0. 



(43) 



(44) 



In this case there are two solutions - again there is one where the propagation direction is parallel to the polarization 
vector — but also another where the propagation direction is orthogonal to the polarization vector (l^v^ — 0). This 
additional solution corresponds to a transversely polarized (vector) sound speed. The two solutions are given by 



f3 + 4/i/3 
P + P ' 



P + P' 



so the two sound speeds are related by 



dP 
dp 



3 Cv 



(45) 



(46) 



For an equation of state where P = wp it can be seen from (|46p that if pj p is sufficiently large then c 2 > even if w is 
negative. This stabilizing property of the shear modulus initially motivated the use of elastic fluids in the framework 
of dark energy models, where w < — 1/3 is required to achieve the observed acceleration [lii , l20l |. The relationship 
between the sound speeds and equation of state is shown in Fig. [21 where we plot lines of constant p/p and c 2 in 
the (w, Cg) plane. The constraint that p/p > requires that w > —1, while the constraint that < c 2 < 1 restricts 
the allowed value of c 2 for a given w. We have assumed that the intrinsic properties of the medium are fixed and 
have ignored time variations in quantities such as w and p/p. Relaxing this assumption could lead to a well-defined 
phenomenological model for time- varying dark energy. However, such a model would have to respect the stability 
conditions discussed above. 
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FIG. 2: Allowed parameter space of elastic fluid models in the (w, cf) plane. The left panel shows lines of constant p/p, with 
the constraint that p/p > disallowing regions with w < — 1 and w > cf . The right panel shows lines of constant cj, with the 
constraint that < (\, < 1 disallowing an additional region with w > — 4/3. 



D. Harmonic Decomposition of Perturbations 

Constant time hypersurfaces are homogeneous and isotropic, so it is natural to decompose spatial tensor fields on 
these hypersurfaces in terms of the irreducible representations of the rotation group SO (3). A spatial tensor field can 
be decomposed in terms of eigenf unctions of the Laplacian [38|, [H, H(| 

V «QW T = -k 2 Q s ' v > T , (47) 

where the S,V,T index represents irreducible scalar, vector and tensorial quantities and '\' represents covariant 
differentiation with respect to the three- metric r/ij. These eigenfunctions form a complete basis in which to expand 
the tensor field. In flat space, for example, Fourier plane waves provide a local orthonormal basis. Scalars can be 
constructed from longitudinal type vectors and tensors via [4l[ 

Qf = -k-'Qfi , Qfj = k~ 2 Q% + ± Vij Q s , (48) 
and similarly vectors can be constructed from solenoidal type tensors by 

QYj = -(kr'Qm . (49) 

subject to Q^ 1 — Q^ 1 — Qf 1 — 0. A vector field can, therefore, be decomposed as 

Si = Z 8 Qf + t v QY, (50) 

and the general decomposition of a symmetric tensor field is 

H l0 = H[Q s Vij + H*Q s t3 + H V QY 3 + H T Q% . (51) 

The six degrees of freedom of the tensor field are represented by two scalar parts along with vector and tensor parts 
each with two degrees of freedom. 
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E. Einstein Equations and Conserved Energy-Momentum 

The perturbed Einstein equations subject to the synchronous gauge conditions are 

-m 2 



a Cx o 

2a 2 G° t 
2a 2 G\ 

a 2 G l , 



— dih — djh 3 i , 
= djh ij -d v h, 



(m - n 2 



(52a) 

(52b) 
(52c) 

(52d) 



djdih'k 



dkdjh) 



The metric perturbation is parameterized by h%j = 2Hij, where Hij is decomposed according to (|51|) . In order to 
compare the elastic fluid with other fluid based models, we use the parameterization of the the energy-momentum 
tensor used in ref. [4lj |. 



T°o = -(fi + Sp), 
T°i = (p + P)v t , 

t\ = -(p+py, 

1*i = (P + 6P)6 i j 



pit 



(53a) 
(53b) 
(53c) 
(53d) 



where vi is the velocity perturbation from the flow and the anisotropic stress, IPj, is symmetric and traceless. These 
quantities are also decomposed according to ([50)) and (f5"Tj) . The scalar- vector-tensor (SVT) split of the perturbed 
Einstein equations then gives the following constraint and evolution equations. 



Constraint: 



Hh-2k 2 r] = 8irGa 2 5p, 

fcj) = 47rGa 2 (p + P)v s , 
-16irGa 2 (p + P)v v 



kH v 



(54a) 
(54b) 
(54c) 



Evolution: 



h + 2Hh-2k 2 r 1 = -2AirGa 2 8P, 
h + 6t) + 2H(h + 6rj) - 2k 2 r] = -16T:Ga 2 PU s . 

H V + 2HH V = 8irGa 2 PIl v , 
H T + 2HH T + k 2 H T = 87rGa 2 Pn T , 



(55a) 
(55b) 
(55c) 
(55d) 



where h = 6-fff and r) = - (H[ + are the metric variables defined in ref. [HJ ■ The equations of motion for the 

scalar and vector components of the displacement vector £ z of the isotropic elastic fluid given by (|35p are then 



e+ i 



3 f) 



3 < k 2 c 2 



2k 
dP 

) 

dp 



2„2 



k z c 



-(H 



0. 



(56a) 
(56b) 



where the subscript / denotes the metric induced on the medium at the time of formation. Assuming an equation of 
state P = wp then the scalar components of the perturbed energy-momentum tensor (|34al I34bl I34c|l are, according 
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to the decomposition (|53al I53bl I53cl I53d[) . 



Sp 



-p(l + w) 



k f + -(h-hi) 



v 
SP 



s ^ 



n 



s 



3(77 - 77/) 



(57a) 
(57b) 
(57c) 

(57d) 



1 + w 



Note that SP/Sp = ^E. Eqn. (|56a[) can then be used in conjunction with (|57al I57bl I57cl l5Tc|) to give a closed set of 
scalar equations 



(1 + w) ( few 5 + -/l 



v s = —H 1 - 3— )v+— kS fcn s . 

\ dp J dp 1 + w 3 1 + w 

The anisotropic stress is zero in the perfect fluid limit (p = 0). Similarly, the vector sources are 

v 

24(1 + w- 1 )(kt v + HY-H v ), 



V t-V 



and the anisotropic tensor source is 



n T = 24(i + w - 1 )( J ffJ-^ 



(58a) 
(58b) 



(59a) 
(59b) 



(60) 



The vector or tensor contributions due to the elastic medium are absent in the perfect fluid limit. Eqn. (|56b[) can 
then be used in conjunction with (|59ai I59b|) to give the vector equation of motion 



- H{ 1-3^ l r 

dp, 



1 w 



2 1 + w 



(61) 



F. Comparison to Generalized Fluid Systems 



The equations of motion derived in the preceding sections are based on a very specific set of assumptions. Given 
that the understanding the evolution of dark energy perturbations is crucial to making precise predictions for the 
observed power spectra, a variety of phenomenological models have been discussed [2II I27I. l29l. l30l. l3lj] . In this section 
we attempt to make some contact between our work and these models. 

Energy conservation of a generalized fluid energy- momentum tensor gives the scalar equations (29l . |42| : 



1 + w 



kv s + i/i 



w 

3H-—T., 

1 + w 



dP\ a dP 1 

N{{ - 3 lp-) v ' 



dp 1 + w 



-k5 



-kT kll s . 



1 + w 



2,1 + w 



where Y is the entropy contribution which is given by 



„r = l"u 

dp dp , 



(62a) 
(62b) 

(63) 



In general, P = wp does not imply SP — w Sp due to temporal or spatial variations in w which correspond to entropy 
perturbations. It transpires that both entropy and anisotropic stress can play a role in stabilizing perturbations which 
is crucial to the viability of any model. 
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In an elastic medium, we have already pointed out that Sp/SP — which implies that T = 0, that is, the medium 
is adiabatic. However, the non-zero rigidity leads to anisotropic stress and it is that which stabilizes perturbations as 
shown by (|45|) . 

A number of authors [29L l30l . l3ll ] have suggested phenomenological approaches which involve the inclusion of 
anisotropic stress. Based on various arguments, they suggested making H s dynamical and constructed a phenomeno- 
logical equation of motion for its evolution 



IT- 



1 



s 



4c 2 



w 



kv" 



3/y 



where T is some decay timescale (which they suggest should be (37Y)" 1 ) and c; 
physical origin is attributed to viscosity within the fluid. If we differentiate 
obtain 



2(i 



kv" 



1 



-h + 3i) 



(64) 

is some arbitrary coefficient, whose 
and substitute in (|58ap . then we 



(65) 



which has a similar form to if T = 00 and c 2 is = fj,/ (2p) = c 2 (l+w)/2. The construction of these phenomenological 
models was somewhat ad-hoc, nonetheless they are very close to those described here. However, we see that the 
physical mechanism which they describe is not viscosity, but rigidity. It might be possible to construct models with 
finite values of T by modifying our approach to include some kind of dissipation. 

We note that the system of equations which we have derived is very similar to that which might come from a 
truncated Boltzmann hierarchy for a black-body such as the neutrinos (or photons with no Thomson scattering 
terms). In particular, if we take w = 1/3 and p/p = 4/15, then the elastic model gives 



-- (kv s + -h 
3 ^ 2 

k (l S -lll s 
W 6 



- 1 



h 

2 



37) 



(66a) 
(66b) 
(66c) 



which is exactly that given by a Boltzmann hierarchy with a third moment set to zero [4l| . 

A second class of models seek to model dark energy fluids which are non-adiabatic and could come from scalar field 
models. For a minimally coupled scalar field with the standard kinetic term 

Sp = <j)8(t) + ^-8(j> , 5P = <j>5(j) - ^-5<j) , (67) 



and hence 5P/8p clearly depends on the frame of the dark energy defined by Scj) and 6cj>. One can choose to work in 
the rest frame of the dark energy defined by (p + P)v l = k l k54>/(t> = 0, where Cg = 8P/5p = 1. More general models 
with non- minimal couplings or exotic kinetic terms could have ^ 1. 

It was suggested in refs. [2^, [2l} to use a system of equations described by (|62ai I62b[) with no anisotropic stress 
but with r ^ 0. In order to take into account the possibility of different frames, the sound speed was deemed to be 
defined in the rest frame of the dark energy and the equations of motion were modified to apply in an arbitrary frame, 
that is, one makes the transformation to 

5 rcst = 5 + 3H(l + w)(v s -B)/k. (68) 

Here B is the space-time component of the metric perturbation and is zero in both the synchronous and Newtonian 
gauges. This leads to equations of motion 



kv" 



1 



1 



k 2 

kc L c 

w 



5 = -(1 + 

i, s = -H{l-3c 2 s )v s 
which can be computed from (|62b[) using 

wT cS = {ci-w)[5 + 3W(1 + 



m(ci - w)s . 



When c 2 — w these are clearly the same as for the elastic medium, but they are very different in the limit c 2 
we shall see in the subsequent discussion. 



(69a) 



(69b) 



(70) 
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Symmetry 


Number of Non-Zero Shear Moduli 


Triclinic 


20 


Monoclinic 


12 


Orthorhombic 


8 


Tetragonal (C4, S4, C41J 


6 


Tetragonal (C 4v , D2d, D4, D 4h ) 


5 


Rhombohedral (C3,Se) 


6 


Rhombohedral (C3 V ,D3,D 3d ) 


5 


Hexagonal 


4 


Cubic 


2 


Isotropic 


1 



TABLE I: Number of non-zero shear moduli for each of the Bravais lattice symmetry groups [431 ]. For the tetragonal and 
rhombohedral symmetries the number of shear moduli depend on the particular class of symmetry, which can be identified 
using the Schonflies notation. There are 7 shear moduli for tetragonal symmetry with a 4-fold rotation axis (C4), for example, 
but this reduces to 6 with the additional of a mirror plane parallel to the axis of rotation (C4 V ). 



G. Anisotropic Generalizations 



In the previous sections we have described how one can construct a general perturbed energy-momentum tensor 
under a set of simple assumptions. The fact that the general isotropic tensor of rank four has a limited number of 
degrees of freedom has allowed us to construct the most general set of equations of motion for an isotropic medium 
defined by its density, pressure and rigidity. More generally we have shown that specification of all of the components 
of the pressure tensor, P^ v , and elasticity tensor E^f" 7 is sufficient to describe the perturbed energy momentum 
tensor using (|12p . In this section we will discuss aspects of a generalization which allows for the perturbations to be 
anisotropic and consider the case of cubic symmetry which is the simplest anisotropic possibility. 

In general, the elasticity tensor has a total of 21 independent components [43| . One of these, the bulk modulus, is 
specified by the pressure and the other 20 are shear moduli, which can provide an anisotropic response to perturbations 
depending on the particular symmetry of the system. Fortunately the classification of these shear moduli has been 
studied in the context of classical elasticity theory and it is known that there 14 different types, known as the Bravais 
lattices. In table U we list the number of non-zero shear moduli for each symmetry group. 

An additional requirement that we will impose here is that the pressure tensor is isotropic, P^ v = Pj^ , so that 
the unperturbed spacetime has the standard FRW metric, and the anisotropic response is only present at linearized 
order. This will probably restrict the range of possibilities allowed, but for sure there is at least one possibility, that 
of cubic symmetry, which is compatible with this. 

The case of cubic response to perturbations was considered in ref. [44| ]. There are now two shear moduli, /xl and 
/it, which specify the non-zero components of the elasticity tensor 



J^XXXX 


= jr;yyyy 


J^ZZZZ 


= [3 + P + 


fixxyy 


= jr;yyzz 


J^ZZXX 


= 0-P- 




j^xzxz 


— J^xyxy 


= P + Ht 



(71) 



The energy-momentum sources are modified from those in (|Mal l34bl IMc]) to 



5T ° = (p + P) (di? + ~h\ , (72a) 

ST o = -(P + P)C, (72b) 
2 



ST] = -8){(i - ~ Mi ) ( 8 k e + - fe(2%f J + h)) - Afj, S l j . (72c) 
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where Afi — /it — Ml quantifies the degree of anisotropy. The tensor S i j represents the cubic source term and is given 
by 



t 

2d {y ^ 



2d {y ^ 
hi 2d {y ^ 

\2d {z ^+h% 2d (z ^+h z y 

are modified to 



(73) 



Furthermore, the evolution equations 

(p + P)(C + HC) - 3j3HC - (3{d l d 3 e + d l h/2) - t i : (<)•<) + d l d^ 3 /3 + 9 s - &h/Z) = A/iF 1 



(74) 



where the cubic source term F l is given by 



F' 



{d y dy + d z d z )£ x + d x {d y £y + d z £, z 



d v h x 



d z h x z \ 



(d x d x + d z d z )£? + dv(d x £ x + d z £ z ) + d x W x + d z W ' z 
\ {d x d x + d v dy)i z + d z (d x (, x + d y £ y ) + d x h z x + d y h z y J 



(75) 



One of the consequences of the introduction of two shear moduli is that the sound speeds become dependent on the 
direction of the wave- vector (45[. When one performs the same SVT decomposition as for the isotropic case, the 
different SVT sectors couple and scalar modes can excite vorticity and gravitational waves [44[ . It is is possible that 
this kind of anisotropic behaviour is responsible for anomalies seen the spectrum of the CMB on very large-scales; 
this is presently under investigation. 



III. ANALYTIC SOLUTIONS 



In this section, we identify the regular perturbation modes which can arise in the early universe in the presence of 
an isotropic elastic medium. This task has been performed in the standard scalar case in ref. [46j]. Observations show 
that the primordial perturbation was most likely dominated by an adiabatic scalar mode generated from fluctuations 
in the metric. However, the most general primordial perturbation can also contain scalar isocurvature modes [46j 
generated from variations in the abundance ratios of different particle s pec ies. In the vector sector, a regular mode 
can also be sourced from non-zero initial photon and neutrino vorticity [47] and primordial magnetic fields can also 
lead to regular modes [48j . 

We consider a universe consisting of cold dark matter (c), baryons (b), neutrinos (is), photons (7) and an elastic fluid 
component (e). We identify specific modes where the elastic medium has an equation of state with w = 0, —1/3, —2/3. 
To identify the regular modes in the early universe we assume that the photons and baryons are tightly coupled due 
to the large Thomson scattering term or- One can then obtain exact equations for the evolution of their velocity. An 
expansion in opacity k^ 1 = an e aT which is valid for max (kn c ,HK c ) <C 1 gives 



■ v 



1 



(1 



R) 
R) 



RHv s ; h 
RHv^ h 



0, 



(76a) 
(76b) 



where R = 3pb/(4p 7 ). In the following we define lu k = CI x Hq, R v — f2„/fi r and i? 7 = f2 7 /fi r . The primordial modes 
are then given by a series expansion of equations in Section [II El in terms of the conformal time r. 



A. Scalar Modes 



In the scalar sector the most general primordial perturbation is specified by a mixture of adiabatic and isocurvature 
modes. Adiabatic coupling between various components in the universe requires that the relative perturbation between 
species, given by 
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vanishes. If the various density perturbations compensate in such a way that the initial curvature perturbation 
(rf) is zero on super-horizon scales then these are termed isocurvature initial conditions. The scalar adiabatic and 
isocurvature modes (for initial perturbations in fluids other than the elastic fluid) are listed in table [Til 

The isocurvature modes which can exist due to an elastic fluid are listed in table HTT1 These modes are interesting as 
they correspond to both initial non-zero density fluctuations and non-zero anisotropic stress, resulting from fluctuations 
of the scalar component of the wordline displacement vector £ . Physically, we expect these modes to arise due to a 
perturbed state of the elastic medium relative to the background at formation, that is, the medium is not formed in 
an equilibrium state. 

Due to our definition of anisotropic stress, the stress term for an elastic fluid, defined in (|57ep . will diverge when 
w = 0. The equations of motion do not suffer any such divergence, as this arises only due to our use of the standard 
decomposition of the energy- momentum tensor in (|53d[) . As such, we list the non-zero initial term for the quantity 
w nf for the elastic anisotropic stress in table IIIII 



B. Vector Modes 

The regular vector mode has non-zero initial photon vorticity, having equal and opposite neutrino vorticity [47| . 
We have computed the initial conditions for this mode in the presence of a elastic fluid, which are shown in table HVl 
neglecting the small contributions from the photon anisotropic stress. 

A regular mode also exists due to a non-zero initial vectorial component of the displacement vector £ v in the elastic 
fluid, which also results in non-zero anisotropic stress (as in the scalar mode) and this mode is presented in table [V] 
In the elastic vector isocurvature mode the photon and neutrino vorticity are zero as they have no source term in the 
tight coupling limit. In both the regular vector mode and elastic isocurvature mode, the elastic fluid decouples from 
the equations of motion in the perfect fluid limit. 

The vector isocurvature solution for aa)=-l/3 elastic component can be found analytically in the radiation era 
due to the absence of photon and neutrino vorticity. This solution is given by 



4w„ 



k 2 + Auj c 



sin 



c v (fc 2 + 4w e )' 



c v (k 2 + 4lo c )- 



4ov 



k 2 + Au> c 



(78) 



It can be shown that this solution gives that presented in table [V] if (J — k 1 when expanded as a power series in r. 



C. Tensor Modes 



There are no tensor modes other than the standard adiabatic mode. The tensor solution can be found analytically 
for a w — —1/3 elastic component in the radiation era assuming that the anisotropic stress of photons and neutrinos 
is zero. Using (|55dj) and (|60| then 



H T = Hj 



(P + 4cJw c )5t 



k 



(fc 2 +4c> e )i 



k 2 + AcIl 



(79) 



This reverts to the standard solution H T — Hfjo(kT) when u c = or c 2 — 0. 



IV. COSMOLOGICAL SIGNATURES 

The observable CMB and matter power spectra were computed by modifying the CAMB software [T(3| to include 
an elastic fluid. We consider the elastic fluid in two scenarios - in one instance it acts as the dark energy component 
in an otherwise standard cosmology. The shear modulus [i stabilizes perturbations when w < 0. We also consider in 
models with Q\ =/= and the elastic fluid as a pressureless component with w = 0, but with a non-zero shear modulus. 
The shear modulus introduces a clustering scale, the Jeans length Aj ~ c s to, where tq is the conformal time today, as 
c 2 = 4/i/(3p) when w = 0. This shares some similarities with a hot (or warm) dark matter component as power will 
be suppressed on small scales, although the shear modulus also allows for vector perturbations in the medium. 
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TABLE II: Regular adiabatic and isocurvature scalar modes. The adiabatic mode results from an initial curvature perturbation (r/) and the isocurvature modes for 
CDM, baryons and neutrinos result from initial fluctuations in density, such that the curvature perturbation vanishes. There also exists a neutrino velocity mode, 
where the neutrino and photon-baryon fluid have a spatially varying relative velocity, balanced so that the net momentum is zero. We list perturbation growth in an 
elastic fluid for these modes. 
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TABLE III: Isocurvature scalar modes due to an elastic fluid having an equation of state with w = 0, — 1/3 and —2/3. The 
isocurvature modes arise from spatially varying fluctuations in the elastic fluid resulting in non-zero initial density fluctuations 
and anisotropic stress. We list the initial perturbation for the parameter combination wllf , due to the definition of anisotropic 
stress in (I53d[) . which is regular even if w = 0. 
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TABLE IV: Regular vector mode with non-zero initial photon vorticity, having equal and opposite neutrino vorticity. HY is 



the first order term in the expansion for the metric perturbation H 
mode. 



We list perturbation growth in an elastic fluid for this 
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TABLE V: Isocurvature vector modes due to an elastic fluid having an equation of state with w = 0, —1/3 and —2/3. This 
mode is the analogue of the scalar isocurvature mode, with a non-zero initial vectorial component of the displacement vector 
resulting in non-zero anisotropic stress. 
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A. Scalar Sector 



The scalar C/s are given by 

Ct = 4tt / d(logk)V s (k)\A e (k,T )\ 2 , (80) 

where Ag(k, r ) is the associated multipole moment for the photon distribution and P s (k) is the initial power spectrum, 
parameterized by V s (k) = Ajfc™ 3-1 , where A s is the initial scalar amplitude and n s is the scalar spectral index. Note 
that V(k) = k 3 P(k)/(2ir 2 ), where P(k) is the matter power spectrum. 



1. Adiabatic Mode 



The adiabatic mode is likely to have dominated the primordial fluctuation as a good fit to both CMB and galaxy 
clustering data can be obtained by this mechanism of structure formation and, assuming a flat universe, specifying 
a total of 6 cosmological parameters - Sib, Sl c , ^4 S , n s , h and the optical depth to reionization kr. A dark energy 
component is required in order to reconcile the observed accelerated expansion [1, [f| d, 0, [1] • 

The introduction of a dark energy component effects the CMB anisotropics both by its influence on the background 
expansion rate and its gravitational perturbations. Any component which has the same equation of state w will give an 
identical overall expansion effect. However, an elastic fluid can potentially be distinguished by the different evolution 
of perturbations compared with other dark energy models. If w < then small-scale anisotropies are unaffected by 
the dark energy component as the perturbations at these scales entered the horizon when the fractional dark energy 
density was negligible. The primary contribution of the dark energy then arises at large angular scales through the 
Integrated Sachs- Wolfe (ISW) effect. The temperature-temperature (TT) contribution of the ISW effect to the scalar 
multipoles is given by 



Cf T =4n J d(\ogk)P s (k) 



d T e- K (<t> + ^)j e (k(T -T)) 



(81) 



where <ft an d *P are the Newtonian potentials [42j and Td cc is the conformal time at decoupling. The presence of 
anisotropic stress alters the evolution of the Newtonian potentials. 

In Fig. [3] we plot the CMB anisotropies for multipoles £ < 25 for the best fitting ACDM model [l|, along with the 
best estimates of the TT WMAP three year data [3, Q . We also plot the anisotropies for a number of elastic fluid 
models, along with the corresponding scalar field model, for various values of w and c 2 . In each case we exploit the 
degeneracy between w and h to rescale the first acoustic peak in order to isolate the ISW effect. It is noticeable in 
both dark energy models that there is a reduction in power at large angular scales for w — —1/3 - this reduction is 
slightly more pronounced for the elastic model due to the anisotropic stress. For the scalar field models, the ISW 
effect is minimised for c 2 = and thereby increases monotonically with c 2 . The elastic model is consistent with this 
behaviour for w = —1/3, but for w = —2/3 c 2 = actually maximizes the ISW effect. 

In Fig. [4] we plot the temperature-polarization (TE) power spectrum for the best fitting ACDM model. The signal 
on large angular scales can be attributed to early reionization. The contribution to the large scale TE power for the 
polarization generated at reionization is given by 



Cj E = 4tt / d{\ogk)P s (k) 



dre K {(j> + i>)ji{k{T -T)) 



(82) 



1"0 O 

-dr — (g(k 2 U + fl) + 2gtl + gU) Je (k(T - r)) 



where g = k e~ K is the visibility function, T re is the conformal time at reionization and the polarization source is given in 
terms of the photon anisotropic stress and the zeroth and second order polarization moments by II = n^/3+Apo+Ap 2 . 
The photon anisotropic stress is the dominant term coming from the free streaming of the monopole at recombination. 
We also plot the TE power spectrum for a number of elastic fluid models, along with the corresponding scalar field 
models, for various values of w and c 2 . We find that precision measurement of the TE cross correlation at low £ can 
provide a potentially powerful discriminator of the dark energy model, if non-zero anisotropic stress is present. The 
variation of c 2 for a scalar field model generating no anisotropic stress has little effect on the TE spectrum. The 



clastic fluid, however, shows significant differences for scales of 
of gravitational potentials during reionization. 



< 10 due to anisotropic stress modifying the decay 
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FIG. 3: (Left) A comparison of CMB TT data with the best-fitting WMAP ACDM model at low multipoles. The parameters 
used are fl c h 2 = 0.104, Q b h 2 = 0.0223, h = 0.734, n s = 0.951, A B = 2.02 x 10" 9 and kr = 0.088 1]. (Right set of 2 x 2) A 
comparison between scalar field and elastic fluid dark energy models at low multipole values. The upper panels show scalar 
field models and the lower elastic models. The left and right panels show values of w — —1/3 and w = —2/3 respectively. In 
each case the angular diameter distance to the first acoustic peak has been rescaled to coincide with the best-fit ACDM model 
to isolate the ISW effect. The appropriate change to the Hubble parameter is h = {0.50, 0.60} for w — { — 1/3, —2/3}. In each 
diagram a dotted line shows the ACDM model. The solid curves (bottom to top apart from the w = —2/3 elastic fluid model) 
are c 2 = 0, 0.2, 0.4, 0.6, 0.8, 1, where c 2 < 2/3 in the w — —2/3 case due to the stability condition (|45|l . 



The possibility of distinguishing two competing models whose differences are only significant on large angular scales 
is limited by the effect of cosmic variance. The probability of distinguishing model A from B, assuming that A is 
correct, is then given by (20| 



P{{ai m }\B)J, A 



i 



a 



(A) 



c, 



(B) 



In 



C, 



(A) 



c, 



(B) 



(83) 



In Fig. [5] we plot this quantity for A = a scalar field model and B = an elastic fluid model in the (w,c 2 .) parameter 
space, and indicate the region where the models can be distinguished at 1— a . We find that this probability is too 
small in almost all of the region of interest to discriminate between models, based on CMB TT anisotropies. In the 
right panel of Fig. [5] we plot the analogous quantity based on TE anisotropies, assuming Kr ^ 0, which shows that 
accurate measurement at £ < 10 vastly increases the area of parameter space that can be differentiated. If kr = 
one would be unable to distinguish between scalar field models from elastic models using TE. 

The large scale structure of the universe depends on the growth rate of perturbations in the various species. An 
important aspect of dark energy is that it has a large Jeans length, Aj ~ c s tq, at the present epoch so that it does not 
cluster on small scales and contribute to measurements of Q m in galaxy clusters. We have assumed that the matter 
power spectrum is modified to 



P 5 (k) = b 2 \5 T \ 



(84) 



where St = J2i This makes this untested assumption that the biasing of the dark energy component is propor- 
tional to its density. In Fig. [S] we plot the matter power spectrum for w = —2/3 elastic and scalar field dark energy 
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FIG. 4: (Left) A comparison of CMB TE data with the best-fitting WMAP ACDM model at low multipoles. (Right set of 2 
x 2) A comparison between scalar field and elastic fluid dark energy models at low multipole values. Parameter values and 
labelling is the same as in Fig(3] apart from c? a where (top to bottom) (? s = 0, 0.2, 0.4, 0.6, 0.8, 1 and Cg < 2/3 in the w = —2/3 
case due to the stability condition (14511 . 




FIG. 5: The discrimination 'landscape' between scalar field and elastic dark energy based on CMB TT (left panel) and TE 
(right panel) data. The dark grey shaded area shows the region of parameter space excluded by the bound < I for the 
elastic component, and the light grey area the region where the two models cannot be distinguished. It becomes difficult to 
distinguish the models in the region of cosmological interest between — 1 < w < — 1/3. 

models, decomposing the power into components due to the various species. 

There are noticeable differences in the two cases. The presence of a clustering scale in the dark energy (of << 1) 
affects the elastic fluid power spectrum in several ways. At small scales, power increases as the sound speed becomes 
non-relativistic. If the assumption (|84[) is correct then an approximate limit of Cg > 10 -3 would be required in order to 
be compatible with large scale structure (LSS) data. It is also noticeable that power in the CDM/baryonic components 
is reduced as the sound speed of the elastic component becomes non-relativistic. In the scalar field dark energy model, 
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FIG. 6: The clustering effect of elastic (left set of 2 x 2) and scalar field (right set of 2 x 2) models as d? — > 0. The solid curve 
shows the total power spectrum for w — —2/3, while the long and short-dashed curves show the decomposition of the total 
power into components due to the clustering of the dark energy and CDM/baryonic species respectively. The panels show a 
dark energy sound speed of cl = 10~ 2 (top-left), cl = 10~ 3 (top-right), ci? = 10~ 4 (bottom-left) and c? = 10~ 5 (bottom-right). 
The dotted curve shows the ACDM model. The total power spectrum has been arbitrarily normalized at k = 10 -3 /i Mpc -1 . 
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FIG. 7: The ratio of the ACDM matter power spectrum and models where the CDM has a non-zero shear modulus, AECDM, 
with sound speeds of c? s — 10 -6 (dotted line), c? s — 10~° (short-dash) and c 2 a = 10~ 4 (long-dash). 



however, a sound speed as low as = 10~ 5 does not change the total matter power spectrum significantly. We have 
reduced of further and find that = would be compatible with LSS data. 

Finally, we also consider whether the introduction of a non-zero shear modulus in a pressureless w — component 
is compatible with CMB and LSS data. In Fig. [7] we plot the ratio of the ACDM matter power spectrum against 
models where there is elastic rigidity in the CDM, which we denote AECDM. The introduction of a non-zero shear 
modulus introduces a clustering scale in the CDM and power is suppressed on small scales. When of = 10~ 5 , for 
example, power is suppressed by an order of magnitude compared to the ACDM model at k = 0.1/iMpc -1 . For the 
values of < 10 -4 , the CMB anisotropics are affected at less than the 1% level by variations in Cg. 
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FIG. 8: Transfer function of the total density fluctuation 5t (left), along with the scale-invariant (middle) and white noise 
(right) CMB power spectrum for the w = elastic model with isocurvature initial fluctuations. The sound speeds are Cg = 



(solid line), 



10 (dotted), c s = 10 (short-dash) and c s = 10 (long-dash). Transfer functions have been arbitrarily 



normalised to unity at k = 1 x 10 5 /iMpc 



2. Isocurvature Modes 



Although the initial conditions appear to be dominated by the adiabatic mode, a sub-dominant isocurvature con- 
tribution cannot be ruled out. The scale-invariant CDM isocurvature mode suffers from the problem that the CMB 
anisotropy has greatly suppressed power on small scales relative to large scales. The reason for this is that the Sachs- 
Wolfe effect and the initial temperature fluctuation both add to give a much greater contribution to large-scale power. 
Recent work has placed an approximate 10% upper bound on the CDM isocurvature fraction, and similar constraints 
also apply for the other isocurvature modes [49j |. 

The power spectrum of the scalar isocurvature mode is defined by 

(|<5| 2 } = / d(logk)Ts(k). (85) 

The scalar isocurvature modes for an elastic fluid are listed in table We list modes for a pressureless w = 
component with a non-zero shear modulus, which can be compared to the standard CDM isocurvature mode, in 
table |TTJ The effect of the shear modulus, and hence the sound speed of, enters at first order in r for several of the 
variables, such as the metric perturbation 77. We also list isocurvature modes where the elastic fluid has an equation 
of state w = —1/3 and —2/3. 

In Fig. [8] we plot the transfer function T(k) of the total density fluctuation St, along with the CMB anisotropies for 
a scale-invariant (Vs oc k , corresponding to P$(k) oc fc~ 3 ) and white noise power spectrum (Vs oc fc 3 , corresponding 
to P$(k) oc k ) for w = and a primordial power amplitude ratio of Vs/V v ~ 1, where Vrj is the power in the 
curvature perturbation. As in the adiabatic case, density fluctuations are suppressed on small scales when c 2 s 7^ 0. 
Since the curvature fluctuation is zero initially, there is the characteristic phase shift of the CMB peak positions 
typical of isocurvature modes compared to the adiabatic one, with the first angular peak at £ ~ 350. For the values of 
Cg < 10 -4 , we again find that the CMB anisotropies are affected at less than the 1% level. When w < perturbation 
growth is suppressed and T{k) is much flatter, with a fall off in power when c 2 7^ 0. This is because as these modes 
cross through the horizon the dark energy density is small, and so the growth in the curvature fluctuation is supressed. 
These are shown in Fig. [51 where we also plot the CMB anisotropies for a white noise power spectrum when w = — 1/3 
and -2/3. 



B. Vector Sector 



Primordial vector modes constitute vortical perturbations in the early universe. The regular vector mode has been 
discussed in ref. (47[, along with a mode sourced by the anisotropic stress of a primordial magnetic field in ref 48]. 
We have also shown in section fill Bl that a regular solution exists with a non-zero wordline displacement £ v , which 
corresponds to non-zero initial anisotropic stress in the elastic fluid. 
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FIG. 9: Transfer function (left) and white noise CMB power spectrum (right) for w = —1/3 (top) and —2/3 (bottom) with 
isocurvature initial fluctuations. The sound speeds are = (solid line), Cg = 0.005 (dotted), cj? = 0.05 (short-dash) and 
Cg = 0.5 (long-dash). 



1. Regular Mode 

In ref. [47| vector perturbations are expressed in the zero vorticity frame, which coincides with the synchronous 
gauge used here. The Einstein equation was written 

k[a + 2H<t) = -87rGa 2 n y , (86) 

where a is the harmonic coefficient of the shear tensor. The relationship H = — ko~ gives the Einstein evolution 
equation (|55cp . and substituting H\ = — kao into the series solution in section UlI Bl returns the same result as in 
ref. [47l ]. The primordial power spectrum is defined by 

(|a| 2 )= J d(\ogk)V a (k). (87) 

In Fig. [TO] we plot the CMB TT anisotropics for multipoles £ < 25 for a scale invariant spectrum {V a c< fc°) with a 
primordial power amplitude ratio of Va/V^ ~ 10 -3 . 

When the vector sound speed is zero the elastic dark energy does not contribute to the anisotropic stress of the 
Einstein equation (|86[) . We find that the large scale CMB power is then the same for both the ACDM model and 
elastic dark energy models with c 2 , = 0. However, w < requires a non-zero vector sound speed otherwise the scalar 
perturbations will be unstable to collapse, as shown by (|l5]l . In the w = —1/3 case, for example, c 2 > 1/4 is required 
for c 2 > 0. However, the effect of the anisotropic stress on the CMB anisotropics at large angular scales is relativity 
small, as shown in Fig. [TO] 
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FIG. 10: CMB TT power spectrum at low multipoles for the regular vector mode with a scale invariant spectrum and primordial 
vector to scalar power ratio of ~ 10~ 3 . The dashed curve shows both a ACDM model and an elastic fluid with w = —1/3 and 
= (corresponding to cjf = —1/3). The solid curve shows an elastic fluid with w = — 1/3 and = d? = 1. 




FIG. 11: CMB anisotropics for the vector isocurvature mode with a white noise ("Pn oc k 3 ) spectrum, w 
(cj = 5/8). The top curve is the temperature spectrum and below this the B- and E-mode polarization. 



-1/3 and c 2 s = 1/2 



2. Isocurvature Modes 



Isocurvature modes can exist in an clastic fluid due to non-zero anisotropic stress. We quantity these modes by the 
spectrum 



<|n| 2 ) = J d(logk)V n (k). 



(88) 



In Fig.[TT]we plot an example of this mode for w = —1/3, with a white noise power spectrum and a primordial power 
amplitude ratio of Vn/V^ ~ 10~ 3 . The CMB anisotropics are significantly suppressed on small scales relative to 
large scales, and are similar to the scalar isocurvature mode. On small scales, the dark energy density is negligible at 
horizon crossing and the resulting growth in a is small. On larger scales, the dark energy density increases at horizon 
crossing resulting in larger growth in a. 
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FIG. 12: Tensor mode CMB TT power spectrum with r = 0.1 at the pivot scale feo = 0.05Mpc~ . The left panel shows the 
spectrum for a w = — 1/3 elastic fluid and the right panel a w — —2/3 model. The labelling of the curves is the same as in 
Fig. [4] and we observe a decrease in large scale power with increasing (? s . 

C. Tensor Sector 

In many inflationary models a nearly scale invariant spectrum of gravitational waves (tensor modes) is produced. 
These models are parameterized by the tensor spectral index n t and the tensor to scalar ratio r = A t /A s , where A t 
and A s are the primordial amplitude of tensor and scalar fluctuations and the tensor power spectrum is given by 
Vt{k) — A t k nt . The dominant tensor source to the CMB TT anisotropy is 



Cj T =4ir I d{\ogk)V t {k) 



dre- K H 



^ je(k(T - t)) 
(fc(ro-r)) 2 



(£ + 2)! 
(£ 2)! 



(89) 



For a scale invariant spectrum as k — > the power in the photon Boltzmann hierarchy remains in the quadrupole and 
thus the tensor contributions have power enhanced on large angular scales. This is shown in Fig. I12( where we plot 
the TT anisotropy for the ACDM model with n t — and r — 0.1 at the pivot scale fco = 0.05Mpc _1 . The elastic fluid 
contributes a source of anisotropic stress to the Einstein equation (|55d[) , and has the effect of damping the evolution 
of H T and reducing the power on large angular scales. This is another potential signature of elastic dark energy 
models. 



POSSIBLE REALIZATIONS OF ELASTIC DARK ENERGY MODELS 



One possible realization of the elastic fluid dark energy models is a frustrated network of non-Abelian cosmic 
strings [T7L fl8L [Hoi l5ll ] or domain walls [2(| H3] ■ In this section we discuss some basic aspects of these scenarios and 
make attempts to link them with the earlier sections. In order to act as dark energy component, any model needs to 
provide significant negative pressure and have a Jeans length which is comparable to the present horizon, which can 
be achieved by having a relativistic sound speed. 

The idea is that a static lattice of topological defects forms at some point after a phase transition in the Early 
Universe (see ref. [53| for a review of the physics of topological defects and phase transitions). Such a lattice will 
have microscopic, mesoscopic and macroscopic scales. The microscopic and mesoscopic scales refer to the core width 
and the lattice cell size respectively; we will not discuss these specifically here and will just assume that the relevant 
theory allows for the formation of a lattice which is stable to physics on these scales. Most important to the present 
discussion is the existence of some macroscopic scale, L. 

On dimensional grounds the density of static cosmic strings scales as p s t r oc L~ 2 and that for domain walls is given 
by Pdw L . The standard assumption of defect evolution based on simulations of the simplest field theoretic 
models is that a self-similar scaling regime is achieved (see, for example, ref. [54[ and references therein) whereby 
L on t. However, if a static, stable lattice forms, then the appropriate macroscopic length scale will scale with the 
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Mt/p 


Ap/p 


SC 


1/6 


1/3 


1/6 


FCC 


2/9 


1/6 


-1/18 


BCC 


0.109 


0.183 


0.074 


Weaire-Phelan (A15) 


0.168 


0.161 


-0.007 


Bergmann (T) 


0.161 


0.161 


< 10~ 3 


Frauf-Laves (C15) 


0.159 


0.166 


0.007 



TABLE VI: Shear moduli values for the cubic symmetry group with different primitive cells. Ap is defined as pt — Pl, and 
is zero in the isotropic case. The last three cases correspond to tetrahedral close-packing (TCP) structures (see, for example 
ref. HI). 



expansion of the universe and L oc a(t). In this case the appropriate values of w are w — —1/3 for strings and 
w = —2/3 for walls and these have been the values which we have focused on in our earlier discussions. It is possible 
that other values of w might be possible if the equation of state of the walls is not Nambu-Goto. For example, the 
energy-momentum tensor of a string with energy density per unit length U and tension T is 

T^(x) = J daS[x - X(a)} (ueX„X„ - jX'^X'^j , (90) 

where e 2 = X' 2 /(l — X 2 ) and X^ — (t, X(<r)) is the position of the string. Using this one can deduce that 

1 r/ T\ _ T~\ 

(91) 



p. - i*. 



'4 «-! 



where (u 2 ) 1 / 2 is the rms velocity of the strings. Hence, one finds that w = —T/(3U) in the static limit. Since causality 
implies that T/U < 1 we have that w > —1/3, with w — in the non-relativistic limit T << U. Similar arguments 
can be put together to show that w > —2/3 for walls. 

The static defect lattice will behave like an elastic solid since it has rigidity. Recently, it has been shown [55[ that 
under the assumption that the lattice is isotropic p/p = 4/15 for static Dirac-Nambu-Goto strings and walls. If the 
rigidity is comparable to this then the lattice will be stable to macroscopic perturbation modes since c 2 and c 2 are 
both positive and the Jeans length will be a substantial fraction of the horizon since c s , c v ~ 0.1. It is clear that an 
exactly isotropic network would be difficult to form in cosmological phase transition since it will be uncorrelated on 
large scales. However, it might be possible for it to be approximately isotropic allowing for the treatment focused on 
in this paper to be applicable. An alternative is that the lattice has approximate point symmetry so the elasticity 
tensor will also have point symmetry as described in section III Gl The shear moduli for the Bravais lattices with 
cubic symmetry (the primitive lattices: simple cubic (SC), face-centred cubic (FCC) and body-centred cubic (BCC)) 
relevant to the domain wall case were discussed in ref. [45| and it was shown that if w = — 2/3 then the BCC lattice 
is unstable and the SC/FCC lattices have zero modes. It was argued that for the SC lattice this mode corresponds 
to perturbation of infinite extent, that is, one of the faces of the cubic system moving toward the another. A number 
of composite lattices with cubic symmetry which are part of the tetrahedral close-packing (TCP) structures (see, for 
example ref. (56[) are also possible, although none of these appear to be stable if w = —2/3. They could, however, be 
stable if w = —2/3 + e for e small and positive. The values of shear moduli are tabulated in table IVll for known cubic 
lattices. 

In addition to not being totally isotropic at formation (or even having exact point symmetry) a lattice may have 
initial macroscopic fluctuations relative to the equilibrium position which correspond to the isocurvature modes 
discussed in the previous section. Since the formation of the lattice is a causal process then these modes with have 
an initial white noise spectrum (Vi(k) oc fc 3 ). 

Two interesting issues are the scale of symmetry breaking of the phase transition and the likely cell size of the lattice 
at the present day. Both of these rely on us understanding the lattice formation process which is not well-understood. 
A simple assumption would be that the lattice forms instantaneously after the phase transition with some initial cell 
size £ c = atf which is some fixed fraction of the horizon size at the time of formation tf. If rj is the symmetry breaking 
scale which relates to the temperature of formation Tf = brj, then assuming the lattice forms in the radiation era the 
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correlation size at formation is given by 

UU) = °-^N-^, (92) 

where Af is the number of relativistic species at formation. In the case of domain walls, if the wall density at formation 
is p(t{) = cr/ 3 /(; c , and assuming the network is subsequently swept along by the Hubble flow, we can estimate r\ to be 

, = 100kev(^) 1/4 (^f (^ 2 ) 1/4 . (93) 

Since the coefficients representing uncertainty in the wall formation process have small exponents, it seems reasonable 
to estimate r\ w 100 keV if fi,j w /i 2 ~ 1. One can also compute the present day cell size in a similar fashion, which 
turns out to be 

/ 3 \ !/ 4 / \f \ -3/8 

£c(M = 100pc(^) (O dw ft 2 )- 1/4 . (94) 



, fo3 / V 100 . 

Since the smallest wavenumber at which non-linear effects become important in CMB codes such as CAMB is k ~ 
0.2 Mpc -1 , it seems reasonable to assume that the linear response to perturbations in such a defect network can 
be treated in the continuum clastic medium framework. The computation of the cell size is more sensitive to the 
uncertainty parameters, but since causality gives the upper limit of a < 1 the continuum medium description seems 
justified. One can also perform a similar exercise in the case of cosmic strings. If the wall density at formation is 
p(tf) = cq 2 /^ 2 then one finds r\ ~ 2TeV and £ c (^o) ~ 1-0 AU, with similar uncertainties in these parameters as in the 
domain wall case. 

The allowed symmetry breaking scale may be increased further if the defect network forms during inflation. In 
this scenario the network is inflated outside the horizon and re-enters during a later epoch thus diluting the initial 
density. If one assumes that e N — Tr/T{ where Tr is the reheat temperature of inflaton and N the number of e-folds 
remaining when the network leaves the horizon during inflation then the symmetry breaking scale of a domain wall 
network is given by 

and the present day cell size is modified to 

W^io-V^f'V^). (96) 

If TV = 30 and Tr, = 1 TeV then r\ w 10 MeV and £ G (to) ~ 100 Mpc. One should note that quantum fluctuations in the 
defect forming field, tp, should satisfy satisfy Sip ~ H < rj, otherwise no phase transition can occur during inflation. 
If V ~ M 4 this restricts the mass scale of inflation to M < 3 x 10 s GeV. 

If a lattice of domain walls or cosmic strings exists, their local gravitational interactions can give rise to further 
observable effects. The symmetry breaking scale rj is limited observationally by the local generation of density 
induced perturbations. In the case of domain walls, the mass per unit area is a ~ rf and so r\ = 100 keV gives 
it~5x 10~ 8 kgm~ 2 , and for r\ = lOMeV the wall has a ~ 5 x 10~ 2 kgm~ 2 . If there are only several walls in our 
current horizon they induce a density fluctuation Sp/p ~ Gato, with a similar sized fluctuations in temperature of the 
CMB. A wall with r\ — 100 keV would induce direct temperature fluctuations of ST /T ~ 10~ 9 , which is well below the 
observed value of 8T/T ~ 10~ 5 . Using this constraint directly, one find that domain walls with r) > 1 MeV would be 
ruled out {57j . However, this assumption assumes that the walls move relativistically, with only several walls in our 
horizon. In the lattice structures we have envisaged here the walls are static and only move with Hubble flow, and 
the CMB distortion in this case can be much smaller (58j . In the case of strings, the mass per unit length is given by 
p ~ rj 2 and so rj = 1 TeV gives p, ~ 9 x 10 -6 kgm -1 . The CMB temperature distortion is of the order ST /T ~ 8irGp, 
and so low energy strings discussed here would not produce significant fluctuations in temperature. 



VI. DISCUSSION AND CONCLUSIONS 



We have studied the cosmological implications of a perfect elastic fluid in the framework of General Relativity. In 
previous work on this subject, Bucher and Spergel derived the equations of motion by variation of the action assuming 
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that the fluid is isotropic [19(. Here, we take a more general approach using the material representation concept and 
obtain equations of motion in terms of the pressure and elasticity tensors of the fluid. This allows us to parameterize 
the most general description of linearized perturbations of the fluid in a compact and transparent manner. 

Under the assumption that the pressure and elasticity tensors are isotropic we derive the Einstein and energy- 
momentum conservation equations using this approach. When the fluid has non-zero rigidity there is a source of 
anisotropic stress which stabilizes perturbations when w < 0, making these models candidates to describe the dark 
energy. The anisotropic stress also interacts with the vector and tensor sectors as the shear modulus induces transverse 
waves in the fluid. This phenomenon is well known in the laboratory. In a non-relativistic (low pressure) deformed 
medium with non-zero rigidity the temperature has both temporal and spatial variations. If the transfer of heat occurs 
slowly then the oscillatory motions in the deformed body are adiabatic. These motions correspond to longitudinal 
(scalar) and transverse (vector) elastic waves, which are related by ci? > 4c^/3 |43j| . In the relativistic (high pressure) 
treatment we enforce adiabicity by assuming that there are no temporal or spatial variations in w and the relationship 
between the wave propagation modes is then given by = w + 4c^/3. 

We find that the elastic fluid model is similar to generalized dark energy models [lij [3(| • The construction of the 
fluid energy-momentum tensor in these models was phenomenological, and it was argued that anisotropic stress can 
be attributed to viscosity in the fluid. We find that the functional form of anisotropic stress suggested in [29|, [3(| is 
identical to our own if the decay timescale associated with the stress term is infinite. In this case, the viscous sound 
speed can be directly related to the vector sound speed of a fluid with rigidity by c^ is = c^(l + w)/2. 

There appears to be two mechanisms for stabilizing dark energy perturbations if we wish to model the dark energy as 
a fluid. Either the fluid can be adiabatic and have rigidity, or it can be non-adiabatic and have entropy perturbations. 
A microphysical realization of the latter possibility is scalar field dark energy. As the evolution of perturbations is 
different in each model, we have investigated the potential observational differences in the CMB and matter power 
spectrum assuming the initial conditions were generated from curvature fluctuations in the metric. We find that 
the possibilities for distinguishing the two models are low using only CMB TT data. However, the TE power at 
large scales is reduced in the elastic fluid models relative to the ACDM and scalar field models with kr ^ 0. This 
increases the parameter space in which the two models can be differentiated, up to w > —0.8. As w becomes closer 
to -1 the dark energy perturbations become less important, and the two models essentially become the same from an 
observational point of view. 

The matter power spectrum outlines an interesting difference between the two models. An important aspect of the 
dark energy is that it has a large Jeans length so that it does not cluster on small scales. In principle though, it can 
cluster on large scales and we have attempted to model this by including dark energy perturbations in Ps{k). As 
c s — > we find that the elastic dark energy contribution to Ps(k) dominates and c s > 10 -3 would be required to be 
compatible with large scale structure data. In the case of scalar field dark energy the contribution to Pg(k) from dark 
energy increases as c s — > but remains sub-dominant even when = 0. 

We have also shown that both scalar and vector isocurvature modes are allowed in an elastic fluid. These modes 
correspond to initial fluctuations of the fluid relative to the background, and give rise to anisotropic stress on super- 
horizon scales. A generic feature of these modes is that power is suppressed on small relative to large scales. The 
dominant mode in these models would have to be the adiabatic curvature mode in order to fit the data but there 
could be a sub-dominant isocurvature component. 

Our approach has shown that we can parameterize the perturbations of the energy momentum tensor in terms 
of the pressure and elasticity tensors of the fluid. In this work we mainly focused on the isotropic case, but this 
formalism can be easily extended if the elastic fluid has point symmetry. The symmetry properties of the elasticity 
tensor can be classified by the Bravais lattices, which are familiar in crystallography. Recently, there has been evidence 
of large scale anomalies in the CMB, most noticeably the North-South power asymmetry [59| and alignment of low £ 
multipoles HH, [12, HH , which appear to persist in the third year WMAP data 64] . A possible solution is that the 
energy momentum tensor of the dark energy is not rotationally invariant at linearized order. We have considered the 
case where the elasticity tensor has cubic symmetry and find that perturbation growth is dependent on the direction 
as well as the magnitude of the wave vector. This is work is developed further in ref. [44j |. and we are currently 
investigating whether these anisotropic models are compatable with the data. 
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